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Observation of a Structure in pp -^ pp^^ near the tttt Threshold and its Possible 
Interpretation by 77 Radiation from Chiral Loops in the Mesonic a Channel 
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The pp -^ ppj'j reaction has been measured at CELSIUS using the WASA 47r-detector with hy- 
drogen pellet target. At Tp = 1.20 and 1.36 GeV, where most of the statistics has been accumulated, 
the 77 invariant mass spectrum exhibits a narrow structure around the tttt threshold, which possibly 
may be associated with two-photon radiation of tt^tt~ loops in the mesonic a channel. 
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The question whether there exists a low-lying scalar- 
isoscalar resonance in the tttt system, the so-called cr- 
meson, has a long history. Starting from the a parti- 
cle problem in nuclear physics, where such an exchange 
particle (possibly in form of a correlated tttt exchange) 
is needed to accommodate for the scalar-isoscalar attrac- 
tion in the A^A^ interaction, the quest for the <j-meson has 
found renewed interest in QCD as the chiral partner of 
the pion — and in this context possibly even as the Higgs 
particle of the strong interaction Jl]. Also the scalar- 
isoscalar meson sector is under much debate presently, 
since there are more states known meanwhile (including 
also possible glueball candidates) than can be fitted into 
a single multiplet. Hence it has been suggested 2], that 
there are actually two nonets, a higher- lying one of pre- 
dominant qq nature and a lower- lying one, encompassing 
the a meson, of dominant qqqq structure. This picture 
of the a would be also close to the result of chiral per- 
turbation theory, where the a emerges dynamically in tttt 
rescattering j^^fil^- 

Experimentally the a is very hard to sense because 
of its expected huge width due to its fall-apart de- 
cay into TTTT. Besides of tttt phase shifts recent experi- 
ments on J/ijj -^ ujTT^TT~ |y], D^ -^ TT~TT^TT^ Q and 
r -^ h'r7r~7T^7T^ 3 h^^^ added evidence for the existence 



of a with a mass in the range 320-480 MeV and a width of 
roughly the same size. To learn more specifically about 
the nature of the a meson, its decay into two gammas 
would be of great help, since this would directly give evi- 
dence for tight particle-antiparticle correlations influenc- 
ing heavily the annihilation radiation j^, [lO? UH uM • ^^ 
fact, measurements |l3,[li| of double pion production in 
e+e~ ^7*7* -^ tttt have been used to get access to the 
<j ^ 77 decay |^ [13 • At energies near the tttt threshold 
these reactions are governed by pion production through 
the Born terms. In a theoretical analysis of these data a 
possible contribution of 7*7* -^ a ^ tttt has been inves- 
tigated and a value of Fcr^^^/Fcr^TrTr ^ 10~^ has been 
extracted [3, [l3 ^^ ^^e pole of the a. 

We have carried out measurements of the reaction 
pp -^ PPJJ at the CELSIUS ring in the energy range 
of Tp = lib - 1360 MeV using the WASA detector to- 
gether with the pellet hydrogen target system [T^. The 
detector has nearly full angular coverage for the detec- 
tion of charged and uncharged particles. The forward 
detector consists of a straw tracker unit followed by plas- 
tic scintillator quirl and range hodoscopes, whereas the 
central detector comprises in its inner part a thin- walled 
superconding magnet containing a minidrift chamber for 
tracking, and in its outer part a plastic scintillator barrel 



surrounded by an electromagnetic calorimeter consisting 
of 1012 CsI(Na) crystals. In these runs triggers had been 
set to allow also for simultaneous measurements of single 
pion, double pion and — in case of the highest energy — 
also 7^-production. In case of pp -^ ppjj the two protons 
have been detected in the forward detector, whereas the 
gammas have been detected in the central detector. This 
way the full four-momenta have been measured for all 
particles, allowing thus kinematic fits for each event with 
4 over constraints. Also in order to be safe from events, 
where particles have escaped detection, we require the 
total energy of the detected particles (before kinematic 
fit) to be equal to that in the entrance channel within 
100 MeV, i.e. within a difference smaller than the pion 
mass. To suppress background further we in addition ac- 
cept only photons with E^ > 50 MeV for the final event 
samples. 

The by far highest statistics has been accumulated at 
Tp = 1.36 GeV, hence we concentrate in the following 
on this energy. Fig. la shows a scatter plot of the 77 
invariant mass M^^ versus the pp missing mass MMpp for 
the selected events. The photon energy resolution is best 
for small energies, hence the resolution in M^^ is best at 
low energies. The situation is reversed in MMpp^ since 
the protons undergo hadronic interactions in the detector 
increasingly with increasing energy. Careful inspection 
of this plot exhibits already evidence for some enhanced 
rate between the spots for tt^ and 77. In order to clean the 
sample from events, where the protons have undergone 
very large hadronic interactions, we employ the condition 
(M^^ - MMpp) /MMpp < 0.4. The emerging projection 
of the scatterplot onto the M^^ axis is shown in Fig. lb. 
A peak-like structure around the tttt threshold appears 
now clearly visible. Finally a kinematic fit ( without 
imposing conditions for the final meson mass ) is applied 
to the sample, resulting (Fig. Ic) primarily in a reduction 
of the width of the peaks, in particular of the r] peak. 
The resulting full width at half maximum is ~ 22 MeV 
for TT^ and r] peaks. The width of the structure at the tttt 
threshold also is consistent with this value. 

An important check whether the peak-like structure 
could be of instrumental origin may be provided by mea- 
surements under different kinematical conditions. To this 
end we have analyzed in the same way data accumulated 
at lower energies and found evidence for this spike also 
at Tp = 1.2 GeV (Fig. Id) despite the lower statistics 
there. 

This spike turns out to be very stable against cuts. 
E.g., increase of the threshold E^ = 50 MeV to E^ = 100 
MeV has no significant effect on this spike. MC simula- 
tions for TT^ and 77 production providing a good reproduc- 
tion of the respective peaks in M^^ give no evidence for 
a scatter of counts (stemming from these reactions) into 
the region between the two peaks. Also from simulation 
of 7r+7r~ production we do not get any contributions in 
the M^^ spectrum. Solely the MC simulation of tt^tt^ 
production does give some scatter of counts into this re- 
gion stemming fioiRppTr^ir^ -^ pp^yy^^ events, where two 
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FIG. 1: (a): Scatterplot of M^^ versus the MMpp for all se- 
lected events at Tp — 1.36 GeV. The dashed line shows the 
cut (M^-^ - MMpp) /MMpp < 0.4 to exclude events, where 
protons have undergone large hadronic interactions in the de- 
tector, (b) - (c): M^^ spectrum (projection of the scatterplot 
above) before (b) and after (c) kinematic fit. The dotted his- 
togram in (c) shows background expected from misidentified 
(d): The same as (c), but for Tp = 1.2 GeV. 



low-energy 7s have escaped detection in WAS A. In fact, 
taking into account the absolute cross section of ^ 200/ib 
for the TT^TrO production at 1.36 GeV [HIl^ about 15% 
of the observed counts between tt^ and 77 peaks are com- 
patible with misidentified ppTT^ir^ events ( see dotted his- 
togram in fig. Ic ). The most dangerous situation for 
producing an artefact peak at 2m^ is, if two tt^s decay 
in such a similar manner that the clusters produced by 
their 7s just merge pair wise into each other. However, 
such a constellation, which is included in the MC simu- 
lations shown in fig. Ic, dotted line, is too rare to pro- 
duce an enhancement at 2m^. We also have investigated 
this special scenario with real tt^tt^ data by deliberately 
merging clusters. Again this operation did not result in 
a reproduction of the observed structure. 



Since none of these simulated processes is able to ac- 
count for the structure observed near the tttt threshold 
and also detailed and comprehensive tests of detector per- 
formance and event structures have not given any hint for 
an artifact, we are led to consider seriously the possibil- 
ity that the observed structure is real and might be due 
to the process pp -^ ppa — > ppjj, in particular also since 
pp -^ pp7T~^7T~ an d pp -^ ppTT^TT^ rcactious are dominated 
by a production |lj, [T^ . 

The statistical significance of the spike observed at 
the TTTT threshold may be estimated directly by inspec- 
tion of Fig. 1. There we see that the spike contains 
roughly Ns ^ 60 counts above a background of Nb ^ 60 
counts at Tp = 1.36 GeV. The significance in standard 
deviations depends on the assumption about the back- 
ground behaviour. Neglect of the statistical fiuctua- 
tions of the background gives S = Ns/^/Nb ^7.7. A 
more reliable estimate is S = Ns/^Ns + Nb ^ 5.5, as- 
suming a fiuctuating background, which is smooth and 
well- fixed in shape. Finally, assuming the full uncer- 
tainty of a statistically independent background results 
Ns/VNs + 2 * TVb ^4.5. Somewhat lower val- 
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1.2 GeV, where 



ues arise, 
Ns^Nb^ 30. 

We note that recently the a ^ 77 decay has also been 
searched for 16] in ^^C(7r~, tt^tt^). The M^^ spectrum 
shown there exhibits some indication of a gentle enhance- 
ment near 2777,7^ . Though this reaction cannot directly be 
compared to our case, it is worth noting that the quoted 
relative upper limit is of similar order of magnitude as 
our observation. 

We next inspect the efficiency and acceptance cor- 
rected angular distributions. For these we reconstruct 
the TT^, T] and a momenta from the respective 7 pairs. 
The angular distributions (j{Qcm) of these reconstructed 
momenta in the overall center of mass system are plotted 
in Fig. 2. For tt^ we obtain an angular distribution close 
to (1 + 3 cos^ 0cm) consistent with A and Roper reso- 
nance excitations in pp -^ ppir^. For r] we see an s-wave 
behavior consistent with the excitation of the /Sii(1535) 
resonance in pp -^ ppri and also in agreement with previ- 
ous measurements [13, U^ • And for the spike region we 
obtain a fiat angular distribution, too, again consistent 
with s-wave production. 

In order to get absolute cross sections we have for con- 
venience normalized our data for the 77 peak to the well- 
known cross section of 4.9/ib for pp -^ ppr] at Tp = 1.36 
GeV [13 , which by use of the known 7^ ^ 77 branching 
ratio of 39% ^ results in a cross section of 1.9 jjh for 
pp -^ ppr] -^ PPJJ' This normalization then leads to a 
cross section value of 4 mb for pp -^ ppir^ which is in 
reasonable agreement with previous measurements ^^. 
For the structure around the tttt threshold the quotation 
of an absolute cross section depends crucially on the as- 
sumption of background beneath this structure. If we 
just take the narrow spike above the smooth continuum, 
then we arrive at roughly 1 — 2/ib for its cross section. 
The acceptance and efficiency corrected as well as nor- 
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FIG. 2: Angular distributions of the 77 pair momentum in 
the overall cms for 7s belonging to 7r° peak, ry peak and to 
the spike region, respectively. The solid lines show fits with 
cr(GcM) '^ 1 + 3cos^GcM for 7r° and (t{^cm) = constant 
otherwise. 
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FIG. 3: Differential cross section in dependence of M^^. Note 
the change of scale at M^^ — 0.18GeV/c^ by more than a fac- 
tor of 100! The dotted curve shows MC simulations for double 
bremsstrahlung adjusted in height to the data. The dashed 
line gives calculations of the diagrams in Fig. 4, whereas the 
solid line represents their coherent sum. 



malized M^^ spectrum is shown in Fig. 3. 

A process which principally contributes to the en- 
ergy region of interest is double bremsstrahlung. Single 
bremsstrahlung has been calculated 21] at Tp = 1.35 
GeV to give a total cross section of about 100/ib for 
Ej > 100 MeV. For double bremsstrahlung we are not 
aware of any such detailed calculation. However, we 
may obtain a first crude estimate for its cross section 
by just multiplying the single bremsstrahlung result with 




FIG. 4: Graphs for the dynamical formation of s-wave vir- 
tual 7T~^7T~ loops generated in the pp collision process and its 
subsequent annihilation into two photons. 



the fine structure constant a. This estimate presumably 
will only be a lower limit, since there are much more 
different diagrams contributing to double than to single 
bremsstrahlung. Anyway from these considerations we 
may expect the double bremsstrahlung cross section to be 
of the same order of magnitude as given by the observed 
continuum of counts between tt^ and 77 peaks. Though 
double bremsstrahlung certainly cannot explain the nar- 
row spike, it appears to be a good candidate for explain- 
ing the continuum beneath the spike. The dotted lines in 
Fig. 3 shows a MC simulation for this process adjusted in 
height to the data and assuming da/dM^^ ~ ^77^ .^^J- 
Qualitatively a cusp-like energy dependence of the 
process pp -^ ppa -^ PPJJ can be obtained by a 
Breit-Wigner ansatz da/d M^^ ^ (e^ + T'^/4)-^ with 
£ = Mjj — rua- and the ( energy dependent ) total 
width r = F^^ + Ft^t^. Below the tttt threshold we 
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rria- ~ Ft^t^ ^ SSOMeF we this way arrive at a qualitative 
description of the observed structure. 

To discuss the possible origin of the spike on a more 
fundamental level, we consider the graphs in Fig. 4, 
where correlations in the a channel are formed dynami- 
cally by pion loops generated in the pp collision process 
and coupled to the 77 channel. The blob represents the 
generation of virtual pion pairs and may be adjusted in 
scale to the total hadronic pp cross section at Tp = 1.36 



GeV. Evaluation of the loop diagrams with a branch cut 
at the 7r+7r~ threshold |23(, leads to a cusp-like struc- 
ture (dashed lines in Fig. 3) with a strength basically 
of order a^ smaller than the total cross section. Though 
we arrive at the proper order of magnitude for the effect 
under discussion, the predicted shape is still quite dif- 
ferent. However, this amplitude is to interfere with the 
one for double bremsstrahlung. With a proper choice 
for the relative phase between both processes we obtain 
a curve, which fits the data remarkably well (solid lines 
in Fig. 3). Note that this cusp effect due to 77 radi- 
ation induced by recombination of chiral loops is very 
general and should appear principially also in channels 
other than the a channel discussed here. 

Summarizing we observe a narrow spike around the 
TTTT threshold in the M^^ spectrum of the pp -^ PPJJ 
reaction. Its angular distribution is consistent with s- 
wave production. In a first attempt for its interpretation 
we consider a channel pion loops , which are generated 
by the pp collision process and which decay into the 77 
channel. Chiral loop calculations of this process reveal 
indeed a cusp-like behaviour, which by interference with 
the underlying double bremsstrahlung background can 
give a reasonable account of the data. If true then this 
observation would be a manifestation of a correlated tttt 
system in the a channel. The correlated tttt exchange 
between nucleons 2J, [23 would not need to be lifted 
above the tttt threshold to get on mass shell, it would be 
so already below this threshold with respect to the 77 
decay channel. Similar loop effects are expected to show 
up also in other channels. 
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